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Present models of the superconducting-to-normal transition in transition-edge sensors (TESs) 
do not describe the current distribution within a biased TES. This distribution is complicated by 
normal-metal features that are integral to TES design. We present a model with one free parameter 
that describes the evolution of the current distribution with bias. To probe the current distribution 
experimentally, we fabricated TES devices with different current return geometries. Devices where 
the current return geometry mirrors current fiow within the device have sharper transitions, thus 
allowing for a direct test of the current-flow model. Measurements from these devices show that 
current meanders through a TES low in the resistive transition but flows across the normal-metal 
features by 40 percent of the normal-state resistance. Comparison of transition sharpness between 
device designs reveals that self-induced magnetic fields play an important role in determining the 
width of the superconducting transition. 
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The sharp change in resistance of a superconduc- 
tor over a narrow temperature range is both a natu- 
ral temperature reference and an attractive thermome- 
ter. A Transition-Edge Sensor (TES) consists of a 2- 
dimcnsional metal film that is electrically biased into 
the superconducting phase transition, where its temper- 
ature and resistance respond to deposited energy [l[. 
TES thermometers have enabled some of the most sen- 
sitive calorimctric and bolomctric measurements known. 
TES measurements of single X-ray, gamma-ray, and al- 
pha quanta achieve the highest resolving powers of any 
energy-dispersive technique: E/AE « 4000-5000 d-Q. 
Arrays of TES microbolometers are integral to modern 
submillimeter and millimeter- wave astronomy, achieving 
microkelvin sensitivity in maps of the cosmic microwave 
background (for example, see Ref. [HI). Despite the broad 
use and success of these sensors, much remains uncertain 
about their behavior, including their internal current dis- 
tribution and the physics that determines the width of 
the superconducting transition under bias. 

There are several models for describing the TES tran- 
sition. Bennett et al. @ recently extended the two-fiuid 
model of Irwin et al. Q where the TES bias is sep- 
arated into supercurrent and quasiparticle components. 
Two-fluid predictions for transition shape were compared 
to data, but the current fractions were treated as fit- 
ting parameters. Sadleir et al. 0] identified the impor- 
tance of the proximity effect and treat a TES as a weak 
link between superconducting leads. Working within the 
weak-link context, Kozorezov et al. Q model a TES 
as a shunted junction. However, it is unclear whether 
a weak-link approach can predict the transition shape 
under bias, and doubtful whether proximitization from 
the leads has a strong role in the large devices studied 
here. The critical-current variation with field that is the 
clearest indicator of weak-link behavior is absent in our 
devices. All three models lack geometric detail; yet, real- 
world TESs have additional normal-metal features de- 



posited in complex geometries on the transitioning film. 
These features are used to reduce noise and control the 
transition width [lo| , and are likely to infiuence the flow 
of the sensor bias current. 

The lack of geometric detail in existing transition mod- 
els means they cannot treat the effects of self-induced 
magnetic flelds from geometry-dependent bias currents 
or the effects of fleld inhomogeneities across a device. 
More generally, these approaches assume that the resis- 
tive surface is a function only of current and tempera- 
ture and neglect external or self magnetic flelds. In gen- 
eral, a nonuniform field creates a distribution of critical 
temperatures and currents within a device that can be 
expected to broaden the superconducting transition. A 
crucial unresolved question is the origin of the transition 
width in a TES, a fundamental parameter that affects 
application-relevant quantities such as sensor dynamic 
range and speed. 

In this letter, we map the current distribution in a 
widely used TES geometry as a function of bias point 
within the transition. Further, we show that geometry- 
specific self-fields affect transition shape and that the 
transition is significantly broadened by self-fields. 

The devices used in this study are representative of 
TESs used for ultrahigh-resolution calorimetry. They 
consist of proximity-coupled layers of Mo and Cu whose 
thicknesses (100 nm and 200 nm, respectively) result in 
a transition temperature Tc for the whole film of ap- 
proximately 121 mK. The bilayer has planar dimensions 
of 350 fim by 350 /im and is deposited on a freestanding 
silicon-nitride membrane. An additional 500 nm Cu layer 
is deposited on top of the bilayer and patterned by use of 
standard photolithographic techniques into banks along 
the bilayer edges parallel to current flow and into seven 
interdigitated bars perpendicular to current flow. In our 
"standard" device geometry, shown in Fig. [T^, both the 
Mo current supply and return lines run along one edge 
of the TES. 
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FIG. 1. Pictures of TESs with two bias return geometries. 
The base sensors are 350 fj,m by 350 fim MoCu bilayers with 
seven normal-metal bars. Current flows from left to right, (a) 
Standard bias return, (b) Meander bias return with galvan- 
ically isolated superconducting Nb trace deposited over the 
TES. For the meander design, the supply and return lines ap- 
proach the sensor as a microstrip. For the standard design, 
the supply and return approach the sensor side-by-side. The 
bar width wi,ar = 16 fim, the bank width Wbank = 17 /xm, the 
bilayer width between bars wu = 32 /xm, the bar length Lhar 
= 300 /^m, the number of normal bars Ni,ar ~ 7, and the bank 
length Lhank = 350 fim. 
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FIG. 2. Estimated field per amp of bias current across our 
TESs for the standard lead geometry evaluated at the film 
centerline (j/o = 0; evaluating Eq. [T] away from the yo — 
centerline reduces the field per amp by a maximum of ~ 8% at 
the film edges xo = ^- and yo — (a) Meandering current 
pattern, (b) Uniform current pattern. The black arrows show 
the current path. Typical bias currents range from 60 fiA to 
180 fj,A. 



We first calculate the magnitude of the self-fields in our 
standard design and assess whether it is reasonable that 
they broaden the superconducting transition. Consider 
the self-field in a film of thickness h and width W, with 
the X-axis parallel and the y-axis perpendicular to the 
plane of the film. Solutions for the current density and 
field profile based on the London equations are given by 
Rhoderiek and Wilson 
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but a uniform current distri- 
bution provides a similar result to the expected self-field 
gradient (see Ref. |ll|, Fig. 1). For simplicity, we assume 
a uniform current distribution within the current carry- 
ing regions and calculate the perpendicular field compo- 
nent at a point {xq, yo) relative to the center of the film 
using: 



By{xo,yo) 



Ho J 



{x + xo) 



2T:{Wh) [x + xoY + {yo - y) 



;dxdy. 
(1) 

By applying Eq. [T] to our standard device lead geome- 
try, we have calculated the perpendicular field across our 
sensors, including the field generated by the supply and 
return leads running along the device edges. We consider 
two potential current-flow patterns in the TESs. One op- 
tion is for current to meander in a serpentine path be- 
tween the normal bars of the device. The second possibil- 
ity is for current to flow as quasiparticles across the nor- 
mal bars so that the current density is uniform through- 
out the TES (horizontally in Fig. [1]). Fig. [2] shows 
the two-dimensional fleld pattern for both a meandering 
current and uniform current. In the meandering case, we 
model the current as traveling in a 30 /im-wide sheet that 
traverses a serpentine path between the seven bars. In 



the uniform case, we assume the current travels in the 
entire 350 jum-wide width of the device. A meandering 
current generates fleld gradients approximately ten times 
greater than those for a uniform current, and the mean- 
dering field pattern is more complex, with several local 
maxima and minima across the device. Furthermore, the 
self-field of the meander dominates the field produced by 
the leads. However, in the uniform geometry, field from 
the leads dominates over the TES current contribution. 

In our devices, bias currents range from 180 /iA to 
60 /iA when the devices are biased between 6 % and 
60 % of the normal-state resistance Rn , and are close to 
100 yuA at 20% of Rn- Combining these values with the 
calculations of Fig. [2] result in a predicted self-field vari- 
ation of 5 /iT-10 for the meander case and ~ 0.5 /xT 
for uniform current flow. 

The most direct effect of a magnetic field on a sensor 
is to suppress both its transition temperature and crit- 
ical current Ic- To quantify this change, we measured 
the reduction in Tc and Ic as a function of externally 
driven perpendicular magnetic field [12|. By combin- 
ing the measured Tc and Ic suppression with the cal- 
culated self-field, we can estimate the contribution of the 
self-field to the transition width of our sensors. Transi- 
tion sharpness is conveniently parametrized by the par- 
tial logarithmic derivatives of resistance with respect to 
temperature at constant current a/ = (To/ Ro)dR/dT\jg 
and with respect to current at constant temperature 
A = iIo/Ro)dR/dI\To- For a 5/iT field, we measured 
a STc = —0.9 mK and a 6Ic = —55 fiA. As a first approx- 
imation, we imagine that gradients in self-field create a 
spread in Tc and Ic within a device. If we assume a lin- 
ear transition, with dR = Rn over the measured STc and 
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Sic, ai ^ {To/Ro)Rn/STc and /?/ - {lo/ Ro)Rn /Sic- At 
a bias of 20 %, aj ^ 670 and /3/ ~ 10. These values of ai 
and /?/ are similar to measured results discussed below 
and presented elsewhere [l3|, confirming that self- field 
magnitudes are sufficient to broaden the superconduct- 
ing transition [l2j . 

To determine the current distribution in a TES, we fab- 
ricated devices with sensing elements identical to those 
of the standard design but with a different return path 
for the bias current. This bias-return geometry is shown 
in Fig. [T}) and consists of a "meander" return geome- 
try, where a 350nm-thick transparent dielectric layer of 
silicon oxide grown by plasma enhanced chemical vapor 
deposition, and a via, allow a 300 nm-thick by 6 /zm-wide 
Nb trace to return the bias current on top of the TES 
in a serpentine pattern between the normal-metal bars. 
When current flows in a serpentine path between the nor- 
mal bars of the device, the meander return geometry pro- 
duces a field that reduces the perpendicular component 
of the magnetic field generated by the bias current. 

Since this magnetic field suppresses superconductivity 
and broadens the superconducting transition, compari- 
son of the transition sharpness between designs can be 
used to probe the nature of the current distribution. It 
is likely that the current flow in the device varies based 
on the bias point. We can see this by considering two 
operating extremes. When the sensor is low in the re- 
sistive transition, the bilayer between the normal-metal 
bars is nearly perfectly superconducting, and current will 
meander between the bars. Higher in the transition, nor- 
mal resistance will appear in parts of the bilayer. Once 
resistance begins to appear, it becomes energetically fa- 
vorable for current to flow across the bars because of 
the greater length of the meander path compared to the 
uniform path, and because the normal resistance of one 
square of bar material {Ra,bar — Rn,bank ~ 4.77m51/n) 
is less than the normal resistance of one square of bilayer 
{Rom = 11.77mf]/n). 

To pursue this hypothesis, we have measured transi- 
tion sharpness at multiple bias points from 6 %-60 % of 
Rn corresponding to bias currents from 180/iA-60^A, 
respectively [l3| . Measured values of a/ are shown in Fig. 
131 The measured values of /?/ show similar trends and are 
therefore presented in the online supplemental material 
(l^ . Low in the transition, below about 30 %Rn, wc ob- 
serve an increase in both a/ and /?/ for the meander bias 
return when compared to the standard geometry. At the 
lowest bias measured, the average values of both parame- 
ters are increased by up to 50 % in the meander return as 
compared to the standard return. The increase in ai and 
/?/ in the meander geometry demonstrates that current 
does indeed travel a serpentine path between the bars at 

low %R]y. 

To understand the results of Fig. [31 we consider the 
resistance contributions of three parallel current paths 
through a TES: (1) resistance from the normal banks 
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FIG. 3. Measured ai versus bias point for the two lead ge- 
ometries. Measurements were taken on four different sensors 
of each bias-return geometry. Multiple sensors were tested 
to account for device-to-device variation, and multiple mea- 
surements were performed with some sensors. In total, there 
are six measurements of ai at each bias point. Plot markers 
indicate average values, thick vertical lines show the standard 
deviation at each bias point, and thin vertical lines depict the 
full range of values at each bias point. Markers for the mean- 
der geometry taken at the same bias as the standard geometry 
are slightly offset to the left along the x-axis for clarity. 



Rbanks, (2) rcsistancc from a meandering path between 
the normal bars Rmeander, and (3) resistance from uni- 
form current flow (bus flow) between the normal banks 
Rhus- Following the definitions of the various geomet- 
ric features given in Fig. [Tl the resistance in the banks 

is given by Rbanks — Lbank /{2wbank)Rn Mank- Rmeander 

depends on the strength / of superconductivity in the bi- 
layer, where / ranges from to 1, corresponding to fully 
superconducting and fully normal. For the devices of Fig. 

E Rmeander 

{2 + Wbar/Wbi + {Nbar - 1)(1 + Wbar / Wbi + 

Lbar / Wbi)} f Ru.bi- Rbus dcpcnds On the distance in the 
bilayer Aq. over which a normal current that has tra- 
versed a bar returns to the condensate as a supercurrent. 
The quantity Aq» is the well known branch-imbalance 



length [14|. Here, Rbus 

2KQ*RuM/{Lbar + Wbi)} - 
Aq. ) + ( A^har - 1 ) - 2 Aq. ) } . 



~ Nbar{wbarRn,ba-r/ Lbar + 
- fRD,bi/i^bar + Wbi){2 {wbt - 

We take device resistance 
R to be the parallel combination of Rbanks, Rmeander, 
and Rbus', when /=!, the device is normal with a resis- 
tance R = Rn. a schematic illustrating the meander 
and uniform current paths in the device and algebraic 
derivations of Rmeander and Rbus are given in the online 
supplemental material [12 1. 

From the expressions above, we expect a cross-over re- 
sistance Rx that depends on Aq. above (below) which 

Gmearider = R^neander ^^SS (greater) than Gbus = -R^J^ ■ 

Values of Gmeander / Gbus are plotted versus R/Rn for 
varying Aq. in Fig. [H Our measurements of transi- 
tion sharpness for the meander and standard geometries 
show that significant current flows in a meander pattern 
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for R/Rn as high as 0.3. Therefore, Gmeander > Gbus 
low in the transition below 0.3R/Rn- At R/Rm > 0.4, 
there is no measurable difference in the transition sharp- 
ness between the two device geometries, indicating that 
the majority of the current is in the bus mode, hence 
Gmeander ^ Ghus at biases abovc OAR/Rn. Low in the 
transition, we take the constraint that Gmeander /Gbus ^ 
0.25. Similarly, high in the transition we take the con- 
straint that Gmeander /Gbus ^ 0.25. Thcsc Conditions 
are met by the curves in Fig. SI where Ag* falls be- 
tween 3 and 6 fim. We note that while our choice of con- 
straints on Gmeander /Gbus IS subjcctive. Varying these 
constraints only shifts the bounds on Aq» . For Ag. 
= 4.5 /xm (used hereafter), Rx/Rn = 0.21. We make 
several observations. First, measurements of a macro- 
scopic parameter, transition sharpness, provide remark- 
ably tight constraints on the microscopic parameter Ag* . 
Second, these results provide clear proof of the impor- 
tance of quasiparticle transport in TES bilayers under 
realistic working conditions, since Ag* = is inconsis- 
tent with the data. Third, because the meander path 
length is greater than the uniform path length and be- 
cause Ru,bar < Ra,bi^ achieved at a very low value, 
fx = 0.005, indicating that the appearance of almost in- 
finitesimal resistance values in the bilayer will begin the 
transfer of current into the bus mode. 

Having determined Ag*, we can compute the fraction 
of the device current in the meander and bus modes, 
as well as in the banks. A current fraction i is given 
by Gi/{Gbus + Gbanks + Gmeander)- The three current 
fractions are plotted versus R/Rn in the inset of Fig. 
131 At low R/Rn, most current is in the meander, but 
current in the bus mode increases with R/Rn, and the 
two modes are equal at Rx/Rn = 0.21. 

The large difference in transition sharpness between 
device geometries sheds light on an important unresolved 
question, namely the origin of the transition width in 
biased TES devices. Not only can the relative values of 
aj and /3/ be used to locate current flow, the quantitative 
differences in ai and /?/ indicate that self-fields play an 
important role in determining the transition width. Since 
transition width is proportional to l/a/, the transition 
in the meander geometry is roughly 2/3 narrower than in 
the standard geometry. Because the designs differ only 
in the self- field experienced by the TES, as much as 1/3 
of the transition width in the standard device geometry 
at realistic bias points between 0.05-0.2 i?Ar must be due 
to self- field effects. 

The calculated self-fields of Fig. [2] provide additional 
insight into the measurements of Fig. [31 Above Rx, the 
current flow switches from a meandering supercurrent to 
uniform current flow. At high bias, we might expect that 
the standard design would have a sharper transition than 
that of the meander design, as the fields from the device 
current and the meander return no longer cancel. How- 
ever, the current is decreasing, and the total field gradi- 




FIG. 4. Calculated ratio of Gmeander /Gbus versus bias 
point as a function of Aq* . The measurements of transi- 
tion sharpness (Fig. [3] and Ref. Fig. 1) indicate 
that Gmeander ^ Gbus for Was values below 0.3 R/Rn, 
but that Gmeander ^ Gbus above biases of 0.4 R/Rn. 
These conditions constrain Ag* to a value between 3 fim and 
6/im, where Gmeander /Gbus ^ 0.25 for R/Rn below 0.3 and 
Gmeander /Gbus ^ 0.25 for R/Rn above 0.4. Insct. Calculated 
fraction of current in the three current paths versus bias for 
Aq* = 4.5 /xm. 

ent across the sensor also drops dramatically, as shown in 
Fig. I^b, so self-field effects simply become less relevant. 

In summary, our findings provide strong evidence that 
device geometry, particularly normal-metal features, in- 
fluences the flow of current in a TES. For a sensor with 
interdigitated normal-metal bars, we have shown that 
current travels in a serpentine path between the bars at 
low %Rn, but that this is completely finished by 40% 
of i?jv • We have developed a model of the current distri- 
bution as a function of device bias and used this model 
to constrain the quasiparticle branch imbalance length 
Ag. . This bias current produces an inhomogeneous self- 
induced magnetic field across the sensor. These self-fields 
contribute to the superconducting transition width. Low 
in the transition, the ai and /3/ differences between the 
meander and standard current returns indicate that self- 
field effects can account for as much as 1/3 of the ob- 
served transition widths. Because TES calorimeters are 
typically biased between 5% and 20% of Rn, self-field 
effects are most relevant at actual working conditions. 
Our measurements of the transition broadening from self- 
fields are likely a lower bound, since several effects con- 
tribute to imperfect field cancellation by the meander 
return: the 350 nm thickness of the dielectric, the some- 
what arbitrary 6 /im width of the Nb meander, and a 
layout error where the meander return enters and exits 
the TES at the top of the device (as shown in FiglT]) but 
should enter and exit at the bottom. 

These results provide insight into the physics govern- 
ing the superconducting transition in 2-dimcnsional films 
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with additional normal-metal features. Our determina- 
tion of the internal current distribution may inform fu- 
ture models that move beyond a strictly monolithic pic- 
ture of TES devices. While we have shown that self- 
fields account for a significant fraction of the transition 
width at low %Rn, other effects must also play a role 
in determining the transition width. A likely mechanism 
is spatial Tc and Ic variation within the device induced 
by the lateral proximity effect between the bilayer and 
normal-bar regions. Finally, we speculate that the ob- 
served suppression of excess noise by normal-metal fea- 
tures is because they rapidly force current into a 
quasiparticle mode that resembles current flow in con- 
ventional resistors and is less susceptible to fluctuation 
mechanisms associated with the superconducting phase 
transition. 
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